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John  D.  Clayton 
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Aberdeen  Proving  Ground,  MD,  USA  21005-5066 
j  clayton@arl.army.mil 

ABSTRACT:  A  continuum  theory  is  developed  for  modeling  elasticity,  plasticity,  and 
twinning  in  single  crystals  of  arbitrary  anisotropy,  subjected  to  arbitrarily  large 
deformations.  Concepts  from  discrete  lattice  models  are  considered  for  estimating  resistances 
to  inelastic  shear.  Nonlinear  thermoelastic  calculations  are  used  in  conjunction  with 
experimental  shock  physics  observations  to  provide  estimates  for  resistances  to  slip, 
twinning,  and  shear  fracture  on  preferred  planes  in  alumina  single  crystals  (i.e., 
rhomobohedral  corundum  or  sapphire). 

INTRODUCTION:  Large  deformation  theory  is  needed  for  an  accurate  description  of 
finite  shear  strains  associated  with  slip  and  deformation  twinning  and  improper  lattice 
rotations  across  twin  boundaries.  Nonlinear  anisotropic  elasticity  is  also  necessary  for 
accurately  describing  the  response  of  ceramic  single  crystals  to  high  pressure  loading,  for 
example  plate  impact  or  explosively  driven  shock  (Graham  &  Brooks  [1971]). 


PROCEDURES,  RESULTS,  AND  DISCUSSION:  The  deformation  gradient  F  is 
decomposed  multiplicatively  as 

} 7“  -  j7Ea  pIa  f pP  -  fe“  Fa  m 

1  .A  1  .a1  .p1  .A  1  .a1  A  ’  V1/ 


where  FE  accounts  for  thennoelasticity,  F1  for  twinning  shears,  and  Fp  for  dislocation 
glide.  Skew  components  of  the  covariant  derivative  of  Fp  can  be  used  to  construct  a  density 
of  geometrically  necessary  dislocations  associated  with  slip.  In  fact,  Nye’s  [1953]  theory 
followed  from  observations  of  slip  traces  in  transparent  corundum.  A  density  of  partial 
dislocations  associated  with  incompatibility  of  twins  follows  from  the  gradient  of  F1 .  A  net 
Burgers  vector  b  is,  from  integration  over  circuit  C  encircling  reference  area  A  with  normal 
N: 


b °  =  -j)c FaAdXA  =  J  daANAdA  , 


daA=sABC{FIapFp^c+FIal}  cFPI,B)  . 


slip  gradients 


twin  partials 


(2) 


Contributions  of  slip  gradients  to  the  Burgers  vector  have  been  described  elsewhere  (Nye 
[1953],  Clayton  et  al.  [2005]).  The  final  term  on  the  right  in  the  second  of  (2)  results  from 
gradients  in  twin  fractions,  e.g.  interface  dislocations  at  tapered  twin  boundaries. 
Disclination  models  of  twins  (Clayton  et  al.  [2005])  could  complement  such  a  description. 
The  Helmholtz  free  energy  per  unit  relaxed  volume  of  crystal  is  prescribed  as 
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•F  =X-Elf.*‘Ell  + 1  Elf.*s‘*E%E% - /T^g*  ($-$„)- c<? log 


second  order  elasticity 
pE 


third  order  elastic  constants 


thermoelastic  coupling 


V^oy 


+  x¥t 


(3) 


lattice 

defects 


specific  heat 

In  (3),  EK  is  the  Green  elastic  strain,  C  are  elastic  coefficients,  p  are  thermal  stress 
coefficients,  and  c  is  the  specific  heat.  Kinetic  relations  for  slip  and  twinning  are  written 


Y  =Ys 


z 

I  z‘ 


y  i  =  rs 


\Tj\ 


Yt 

A) 

sJ 

gJ 

(4) 


where  y  and  ;/'  are  slip  rates  for  slip  system  i  in  the  parent  crystal  and  for  slip  system  i  in 


twin  fraction  j,  f  J  is  the  rate  of  evolving  volume  fraction  of  twin  j  with  characteristic  shear 
sJ ,  z  is  the  resolved  Eshelby  or  Kirchhoff  stress,  g  are  resistances  that  evolve  with  the 
thermodynamic  state,  and  remaining  terms  are  constants.  The  crystal  structure  and  relevant 
systems  for  slip  and  twinning  in  sapphire  are  shown  in  Fig.  1 .  A  full  close-packed  layer  of  O 
atoms  lies  in  between  each  partially  filled  (0001)  plane  of  A1  atoms  in  the  hexagonal  cell. 


Fig.  1  Unit  cells  of  sapphire,  after  Kronberg  [1957],  and  slip  and  twin  systems. 


Because  of  sapphire’s  brittle  nature,  values  of  g  must  be  inferred  from  experiments  where 
confining  pressures  are  high,  such  as  shock  loading  (Graham  &  Brooks  [1971]),  indentation 
(Tymiak  &  Gerberich  [2007]),  or  must  be  extrapolated  from  high  temperature  data.  Previous 
research  indicates  that  the  Peierls  stress  concept  may  provide  insight  into  to  the  low 
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temperature  slip  resistance  (Farber  et  al.  [1993]);  the  ratio  of  stacking  fault  energy  to  Burgers 
vector  may  approximate  the  twin  nucleation  stress;  and  the  theoretical  shear  fracture  strength 
of  the  crystal  may  be  approached  in  shock  physics  experiments  (Graham  &  Brooks  [1971]). 
Values  of  these  theoretical  quantities  are  found  as 

T^er,s  =(2jU/K)exp(  -2nd  /(Kb)),  r^in=WSF/b,  Tf™c,ure  =  jub  /  (2nd) ,  (5) 

where  //  is  the  rhombohedral  shear  modulus,  K  accounts  for  dislocation  line  geometry  and 
anisotropy,  b  is  the  magnitude  of  the  local  Burgers  vector,  d  is  the  interplanar  spacing,  and 
WSF  is  the  appropriate  stacking  fault  energy.  Anisotropic  nonlinear  elasticity  is  used  here  to 

compute  more  accurate  estimates  for  rc ,  and  hence  g  in  the  limit  of  rate  independent 

inelasticity.  Resolved  shear  stresses  for  basal  and  prismatic  slip,  basal  and  rhombohedral 
twinning,  and  fracture  on  pyramidal  and  rhombohedral  planes  are  calculated  during  uniaxial 
straining  of  single  crystals  of  various  orientations.  Stress  values  computed  at  the  Hugoniot 
elastic  limit  (Graham  &  Brooks  [1971])  provide  the  bounds  in  the  first  data  column  of  Table 
1.  These  bounds  compare  favorably  with  those  from  indentation  (Tymiak  &  Gerberich 
[2007]),  but  differ  substantially  from  order-of-magnitude  estimates  computed  via  (5). 


Table  1:  Bounds  on  low  temperature  shear  strengths  for  slip,  twinning,  and  fracture. 


Mechanism 

rc  [GPa]  present  work 

rc  [GPa]  theory  (1) 

tc  [GPa]  indentation  n 

basal  slip 

2.7  <  tc<  5.0 

9.0 

4.8  <  tc<  23.7 

prism  slip 

2.3  <  rc  <  5.7 

0.3 

2.5  <  rc  <  4.1 

pyramidal  slip 

- 

0.7 

7.0 

basal  twin 

rc  >  4.0 

7.3 

4.8  <rc<  7.3 

rhombohedral  twin 

1.0  <  rc  <  8.1 

3.6 

1.0  <  rc<7.3 

prism  fracture 

rc>5.9 

40.5 

- 

pyramidal  fracture 

rc  >  6.7 

30.6 

- 

rhombohedral  fracture 

^  >  7.9 

34.4 

~ 

(1 ’Equation  (5)  <2)Tymiak  &  Gerberich  [2007] 
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